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ABSTRACT 

We study Kelvin-Helmholtz (KH) instability at the interface of a disc and corona sys- 
tem by doing a linear perturbation analysis. The disc is assumed to be thin, however, 
the corona is considered to be nearly quasispherical because of its high temperature. 
Under these circumstances, the interface is subject to the KH instability for a given 
set of the input parameters. Growth rates of the KH unstable modes are calculated for 
a wide range of the input parameters. We show that for a certain range of the pertur- 
bations, the unstable KH perturbations are growing with time scales comparable to 
the inverse of the angular velocity of the accretion disc (dynamical time scale). Thus, 
KH instability at the interface of a disc-corona may have enough time to affect the 
dynamical structure of its underlying accretion disc by possible exchange of the mass, 
angular momentum or even energy. Our linear analysis shows that KH instability may 
provide a mechanism for such exchanges between a disc and its corona. 
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1 INTRODUCTION 

The standard t h eory of the accretion discs 
l|Shakura fc Sunvaevl 1 19731 ) has been used widely over 
recent three decades. Assuming the disc is radiating like a 
black body object, we can explain the thermal emission of 
some of the accretion discs using the standard theory. How- 
ever, non-thermal emission is also detected from some of 
the accretion discs, namely discs around compact objects or 
even active galactic nuclei (AGN). One can hardly explain 
the non-thermal emission based on the standard theory of 
the accretion discs and for this reason non-standard mod- 
els a re proposed (e.g., iKawaeuchi. Shimura fc Mineshigel 
2001). In fact, X-ray of many accreting binaries and 
AGN imply a hot region around the system and these 
considerations have motivated some authors to propose a 
two-component model for these accreting systems, in which 
a thi n disc is surrounded by a hot and tenuous corona 
(e.g.. lHaardt fc Maraschl 1 199 it iNakamura fc Osakil Il993t 



rotational instability (MRI) is the main driving mechanism 
of th e turbulence in the accretion discs ijBalbus fc Hawlevl 
Il99ll ). Numerical simulations of MRI show the magnetic 
field lines may come out of the disc because of the buoyancy 
force and generate the magnetic loops outside the disc and 
dissipate th e energy, mainly through the reconnection mech- 
anism (e.g.. iMiller fc Stondl200(]| ; iLiu. Mineshige fc Shibatal 
2002). Then, the region around the disc is heated signifi- 
cantly bec ause of the magnetic energy dissipation of the 
loops (e.g.. iHewaerts fc Priest! 1 19891; iField fc Rogers! Il993l : 



iBlackman fc Fieldll2000l : lKuncic fc Bicknellll2004) . However, 
there is also another scenario for the formation of a corona 
around an accretion disc. Some authors believe that the 
coronae are created by the radiative pressure driven winds, 
which produce a hot plasma above an accre tion disc (e.g., 
lOhsuga. Mori. Nakamoto fc Mineshigdl2o"05h . 



iKawaguchi. Shimura fc Mineshigdl200ll ; iMerlonil 12003 ). On 
the other hand, it is now commonly accepted that magneto- 
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Since the first phenomenological model of 

lHaardt fc M araschi (|l99ll ). the model has been ex- 
tended in various directions. However, it is still uncertain 
if the coronae are formed because of the magnetic field 
or radiation driven winds or even a combination of both 
mechanisms. Whatever is the formation mechanism of a 
corona, the existence of a hot gaseous system above the disc 
will modify the structure of the underlying accretion disc 
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because of the mass, energy and even angular momentum 
exchange between the disc and corona. Although the 
true nature of these exchanges is yet to be understood, 
one can consider physical arguments to present models 
for disc-corona systems including possible exchanges. 
The basic work about the mass exchange between the 
disc a n d corona was given by Meyer fc Mever-Hofmeister] 
111994 ). IRozaNska fc Czernvl (|2000h studied a hot advec- 
tive corona in pressure equilibrium with the cold disc. 
In their model, the corona exchanges energy with the 
disc through t he radiative interaction and cond uction. In 
another study, Maciolek-Ni edzwiecki et al.l (|l997l ) presented 
a model for the structure of a hot corona by considering 
the possible effects of thermal conduction. Most of the 
semi-analytical models for disc-corona systems a r e following 
a sim i lar approach (e. g. ; RozaNsk a et al.l 1999; IRozaNska 
19991; iLiu et al.1 |2002| ; iMever-Hofmeister fc Meverl |2006| ; 
Mayer fc Pringld 120071 ). Recently, iMaver fc Pringlef (|2007h 
presented time-dependent simulations of a two-phase 
accretion flow around a black hole. They also considered 
mass and energy exchanges between disc and corona. 

In most of the previous studies, for simplicity, it is as- 
sumed the corona is rotating with a Keplerian rotational 
velocity like the disc. Also, the geometrical shape of the 
corona is not investigated in detail. However, the temper- 
ature of the corona is much larger than the disc because 
the observed high energy cut-off in X-ray spectra usually 
is about 100 keV. Moreover, the corona is not thin com- 
paring to its cold disc and the gradient of the pressure in 
the corona can not be neglected. Therefore, in our analy- 
sis, a corona tends to have a quasi-spherical shape with a 
non-Keplerian rotational velocity because of the gradient of 
pressure. However, in the underlying cold disc the gradient 
of pressure can be neglected, and it rotates with Keplerian 
velocity. In other words, not only there is a significant den- 
sity contrast between the disc and corona, a velocity drift 
should exist between the disc and corona. 

When two fluids are in relative motion on either side of 
a common boundary, the Kelvin-Helmholtz (KH) instability 
can occur. Thus, we lead to a simple question: Is the interface 
of a disc and corona system unstable because of KH insta- 
bility? Our goal is to address this question in this paper. KH 
instability has already been studied in different astrophys- 
ical systems from the interactions between the solar wind 
and planetary magnetos pheres w hich are studied based on 
the KH instability fe.g..jNaganolll979T l to jets and outflows 
(e.g., [Hardee 1997; Rossi 1997h and dusty layer in protoplan- 
etary discs (e.g., Michikoshil 12006) . Possible instabilities at 
the interface between an accretion disc and a magnetosphere 
surrounding t he accreting mass ha ve been studied by many 
authors fe.g.. lLTfc Naravan| [2004). It was shown the mag- 
netospheric interface is gener ally Rayleigh- Taylor unstable 
and may also be KH unstable (|Li fc Naravan|[2004 ) . But the 
KH instability at the interface between a disc and its corona 
has not been studied as far as we know. In the next section, 
general equations of our model are presented. Analysis of the 
linear perturbation analysis is done in section 3. A summary 
of the results are discussed in section 4. 



2 GENERAL FORMULATION 
2.1 A disc and corona system 

We consider a rotating reference frame at a distance R from 
the central object (e.g., a black hole or a neutron star) in 
which the corona is at rest. If we denote the angular ve- 
locity of the corona by fl c , the continuity and momentum 
equations in this frame become 

^+V.(pv) = 0, (1) 
<9v VP 

— + (v. V)v = V* — 2fl c x v — n c x (Q c x R) , (2) 

where p, v, P are density, velocity and pressure, respectively. 
Also, ip is the gravitational potential due to the central ob- 
ject. The coordinate system (x,y, z) is oriented so that with 
e z parallel to f2 c and e x and e y pointing in radial and az- 
imuthal directions, respectively. 

Since the corona is geometrically thick, radial pressure 
gradients leads to a deviation form pure Keplerian rotation 
(Mayer & Pringle 2007). We can write the x component of 
the momentum equation @ as 



vL, 1 dP r 

H p c oR 



0. 



(3) 



where Q.k = y GM/R 3 is the Keplerian angular velocity. 
Thus, 



n e _ / dPJdR 

~ V + n 2 K R Pc ' 



V 7 ! - {H c /RYxp ' 



(4) 



(5) 



where \p = 9 (log P c )/<9(log R) and H c is the scaleheight of 
the corona (Mayer & Pringle 2007). 

Now, we can calculate the rotational velocity of the disc 
in our rotating reference frame. Since the disc is cold, we can 
neglect the gradient of the pressure in the x component of 
equation of motion ([2]). Thus, 

2n c v 0y + {n 2 c ~n 2 K )R = o, (6) 

and using equation Q, the above equation can be written 
as 

PcXP 



V y = 



2Q c Rp c 



(7) 



Q C H C , we get 



(8) 



Since P c ~ p c c StC and c s 

y 0y = -7^c(^t) 2 xp- 

Substituting f2 c from equation ((5]l into equation ([8]), we ob- 
tain 

9 



V y 



-VK, 



(9) 



2V1TB 

where O = —(H c / R) 2 \p ls a nondimensional parameter 
which depends on the physical structure of the corona. This 
parameter shows deviation of the rotation of the corona from 
Keplerian angular velocity. Obviously, the parameter O de- 
pends on the physical properties of the corona which are ob- 
tained via solving dynamical equations of the corona. Our 
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Figure 1. Schematic of the discrete three layer model for a disc 
and corona system. Generally, thickness of the corona is much 
larger than the thickness of the disc. Although the disc is rotating 
with Keplerian profile, the rotation velocity of the corona is not 
Keplcrian because of the gradient of the pressure in the corona. 
In the reference frame rotating with the corona at O c , the disc 
has velocity V at radius Rq. 



goal is not to study the detailed structure of corona, though 
it is a difficult task. Instead, we consider the parameter Q 
as a free parameter of our model to study possible instabil- 
ities at the disc-corona interface. If we neglect the gradient 
of the pressure of the corona, then B = and the disc and 
corona are both rotating with Keplerian angular velocity. 
Obviously, in this case, the KH instability is not important. 
Because both disc and corona are rotating with the same 
Keplerian angular velocity and there is not drift velocity be- 
tween them. Thus, in our work, we consider non-zero values 
for 6. 



2.2 Linear perturbations 

We assume incompressibility for the analytical calculations 
and all unperturbed physical quantities are considered con- 
stant in both disc and corona. We perturb the density, the 
pressure and the velocity of the system as 



p{z,y,t) = pi{z) exp[i(kyy - ut)], 
P{z,y,t) = P 1 {z)exp[i{k y y - ut)], 
v(z,y,t) =vi(z)exp[i(k v y-ut)]. 



(10) 



After doing mathematical manipulations, the final differen- 
tial equation for the perturbed z component of the velocity 
becomes 



V{p(cj — k y voy)[l — 



(w — kyVoy) 



kyp{u> — k y v y)vi z = gkyVp- 



(11) 



where T> = d/dz. In our rotating reference frame, the corona 
is at rest and the disc has a constant velocity which is given 
by equation ©, i.e. 



v = , e v K . (12) 

Figure [T] shows schematic of the discrete three layer 
model for a disc and corona system. The differential equa- 
tion pip is valid at three layers except at the boundaries 
+Hd and —Hd- Since the density of each layer is assumed to 
be constant, the general solution of the differential equation 
(|1H) becomes 



(*0 



de- qiz 
C 2 e +q2Z + Cae" 
Cie qiz 



for z > Hd 
for \z\ < H d 
for z < —Hd 



where Ci, C2, C3 and C4 are the constants of integration yet 
to be obtained from the boundary conditions appropriate to 
our system. Note that qi and qi are given by the following 
equations, and are assumed to have a positive real part so 
as to render the perturbations bounded at infinity. We have 

qy = ky _, = , (13) 



52 = 



4^2/^2 



^/l-mi/iuj-kyvy 



(14) 



In general, we may have solutions which are either 
midplane-symmetric or midplane-antisymmetric. We obtain 
the dispersion relation for both the cases. 



2.2.1 Even solutions 

If the solutions are even in 2, we have v\ z {z) — Vi z (—z) 
which implies C\ = C4 and C2 = C3. Applying the continu- 
ity of the vertical displacements across the interface z = Hd 
and integrating equation (|11[) over an infinitesimal region 
around the interface and use standard Gaussian pillbox ar- 
guments, we obtain the final dispersion relation as 



A(Lj-k y Vf[l 



4ni 



(uj-k y vy 



it/2 



tanh(q 2 -Hd 



+ ^ (1 _4S^ r /2 = 



(15) 



where A = pd/pc and g = g(H d ) = ti 2 K H d . 
2.2.2 Odd solutions 

In this case, the solutions are midplane-antisymmetric, i.e. 
viz{z) = —vi z (—z). So, we have C\ = — C4 and C2 = — C3 
and finally 



A{u-k y Vy[l 



lV2, 



3th(g 2j ff d ) 



+ cAl-^) 1/2 = <?MA-l). 



3 ANALYSIS 



(16) 



To study the nature of the solutions of our nonalgebraic 
equati ons, one can follow a graphical procedure similar to 
that in IChandrasekliarl |l98ll . p. 498). However, we solve the 
dispersion relation numerically to find unstable modes and 
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Figure 2. Growth rate of the unstable modes versus the 
wavemimber for A = 10 3 , = 0.1 and H d /R = 10 -3 . Five dif- 
ferent regions are marked by numbers. The colours are different 
roots. 



their growth rate. Notice that both the even and odd modes 
have the same instability condition for short wavelengths, 
i.e. kyHd 2> 1. In this case, the dispersion equation (|15|) sim- 
plifies to the same dispersion relation as the two-lay er prob- 
lem with rotation studied by IChandras ekhar (198l|). Thus, 
our analysis is restricted to the long-wavelength regime, i.e. 
kyHd <C 1- In the following we present the dependence of 
the unstable even solutions on the input free parameters of 
the analysis. However, there a point that we should note. If 
A represents the wavelength of the perturbations, then we 
have A < 2nR. Thus, k y H d > Hd/R. Although the plots 
are given for a fixed range of k y Hd, we should note that the 
acceptable range is Hd/R < k y Hd < 1 and it depends on 
the opening angle of the disc. 

Now we can introducing the nondimensional growth 
rate F and wavenumber x as 



T = — , x = k y H d . 



(17) 



Therefore, dispersion equation (|15p for even solution is writ- 
ten as 



p 4 r 4 + p 3 r 3 + p 2 r 2 + p x r + p = o, 

where the coefficients are 

p _ [e 2 A:r 3 -4(e + l)(A-lM^) 2 ] 2 
16(9 + 1)2(^)4 

QKx i [9 2 A:r 2 - 4(6 + 1)(A - l)^) 2 ] 



(18) 



Pi 



2(e + i) 3 / 2 (-|i) 3 



36 2 A 2 x 4 



Ps = 



2(e + i)(^) 2 i + e 

29AV 



2(A - l)Ax 2 



P 4 =A 2 x 2 - 1. 



(19) 

(20) 

(21) 
(22) 
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Figure 3. Growth rate of the unstable modes versus the 
wavenumber for different set of the input parameters which are 
shown on each plot. 



There are three input parameters, A, Q and Hd/R- Consid- 
ering complex dependence of the coefficients on the param- 
eters, it is unlikely to give analytical roots and classify the 
stability region in the parameter space in closed mathemat- 
ical forms. However, we solve the above algebraic equation 
numerically for the parameters within the acceptable ranges. 
Roots with non-zero imaginary part are corresponding to the 
unstable modes. But, we accept only those roots with non- 
zero imaginary part which are giving k y Hd less than unity 
and qi and q2 with positive real parts. 
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Figure 4. The same as Figure [3] but for another set of the input 
parameters. 



Figure [2] shows typical profile of the growth rate ver- 
sus wavenumber for A = 10 3 , 6 = 0.1 and H d /R = 1CT 3 . 
The colours are actually different roots. There are different 
regions which are labeled by numbers. In region 1, there is 
only one unstable mode up the wavenumber k y Hd = 0.001. 
After this wavenumber, we can see two unstable modes (re- 
gion 2). Growth rate of one mode is slightly larger than the 
other mode in this region. In fact, the transition between 
the regions 1 and 2 is determined when the coefficient Pa 
becomes zero, i.e. Xt = k y Hd = 1/A. As we will show for 
different values of the input parameters this is the case. In 



each of the regions 4 and 5, there is only one unstable mode. 
Also, there is a stability region between regions 4 and 5 for 
this set of the input parameters. This overall pattern of the 
growth rate versus the wavenumber is seen for a wide range 
of the input parameters as we will see. Interestingly, growth 
rate becomes comparable or even faster than the dynamical 
time-scale for a certain range of the wavelengths. 

In Figure [3] we keep the parameter O fixed at 0.1 but 
the rest of the parameters are varied. We can see the transi- 
tion wavenumber x t shifts to the smaller values as the ratio 
A increase. So, the regions 1 and 2 shift to the left as the 
ratio of the densities of the disc and corona increases. How- 
ever, the location of the region 3 is independent of the ratio 
A, though the growth rate in this region increases with the 
parameter A. On the other hand, while the stability region 
4 vanishes for A = 100, the system becomes more stable in 
this region as the disc becomes more dense in comparison 
to the corona. When the opening angle of the disc increases, 
all the regions shift to the right so that the stability region 
is beyond the acceptable range of the wavelengths of the 
perturbations. 

In Figure [4] the ratio A is kept fixed and the rest of 
the parameters are changed. As the parameter increases, 
not only the curves are shifting to the left, the stability re- 
gion 4 becomes smaller. Figure [5] shows growth rate versus 
wavenumber for a wide range of the input parameters. We 
can see the overall behavior is similar to the previous figures. 
In general, increasing A, Hd/R shift the plots to the right 
and increasing shift to the left. There are some prominent 
features, which are altered as the parameters change. At 
x — 1 /A, there is often a very sharp peak comparable to the 
dynamical time-scale. This special wavelength can be writ- 
ten as A = 2nHd(pc/pd)- There is a range of wavelengths for 
which there are no unstable solutions, and thus the system 
is stable. 



4 DISCUSSION 

We have shown that the interface between a disc and its 
corona is unstable subject to the KH instability and the 
growth rate of the unstable modes was calculated. Different 
time scales are defined for an accretion disc like dynamical 
and viscous time scales. In most of the analytical models of 
the accretion discs, dynamical time scale is less than viscous 
time scale. For some of the input parameters and a range 
of the wavelength of the perturbations, the time scale of 
growing KH unstable modes is comparable to the inverse of 
the angular velocity of the disc (i.e., dynamical time scale). 
It means that KH instability has enough time to affect the 
structure of the disc. Growth rate of the perturbations versus 
wavenumbers shows a very complicated profile depending on 
the input parameters. However, when the wavelength of the 
perturbations is around A = 2nHd(p c / pd) the time scale of 
growth is of the order of the dynamical time scale. 

It is discussed that the cor ona may form because of a 
large scale magnetic field fe.g.. [Blackman fc Pessahll2009l ). 
The exact amount of the energy transfer from the disc to the 
corona by the magnetic flux tubes is yet to be understood 
theoretically, though the re are some est i mates based on the 
physical considerations l|Merlonil 120031 ; iKuncic fc Bicknelll 
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Figure 5. Growth rate of the unstable modes versus the wavenumber for a wide range of the input parameters. Like Figures [3] and [4] 
black, red and green lines are corresponding to A = 10 2 , 10 3 and 10 4 , respectively. Parameter © varies from 0.1 to 3 and the opening 
angle of the disc, i.e. H^/R changes from 0.001 to 0.01 



12004 iBlackman fc Pessahl [20091 '). This possible exchange of 
the energy not only modifies the structure of the corona, 
but may have significant effect on the structure of the disc 
as well (jKhaienabi fc Duffvl 12008). Our analysis proposes 
another mechanism to be considered for the exchange of the 
energy or mass between the disc and corona. Evidently, our 
analysis does not allow to estimate possible momentum and 
energy exchanges between the disc and corona because of 
KH instability. Further numerical simulations are needed to 
address this issue. 
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